Spinless Hartree-Fock model of persistent current in rings with single scatterer: 

Comparison with correlated models 
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Using the self-consistent Hartree-Fock approximation for spinless electrons at zero temperature, 
we study the persistent current of the interacting electron gas in a one-dimensional continuous 
ring containing a single 8 barrier. We calculate the persistent current as a function of the ring 
circumference, magnetic flux threading the ring, barrier strength, etc. We compare our results with 
the results of correlated models like the Luttinger liquid model and the Hubbard model solved by 
means of the renormalization group. A good agreement is found. First, the persistent current decays 
with the increasing ring circumference (L) faster than L _1 and eventually like L _Q_1 , where a > 
depends only on the electron-electron interaction. Second, the persistent current is a sine-shaped 
function of magnetic flux. This sine-like dependence and in particular the universal power law 
j^—a—i jj ave so f ar been believed to arise only in the correlated many-body models. Observation of 
these features within the Hartree-Fock model is a surprising result. 

PACS numbers: 73.23.-b, 73.61.Ey 



I. INTRODUCTION 

Electron gas in a quantum wire is a realistic one- 
dimensional (ID) electron system. By tying the wire ends 
to each other one creates a ID ring. If the ring circum- 
ference is shorter than the electron coherence length, one 
speaks about the mesoscopic ring. An external magnetic 
flux applied through the opening of such ring gives rise 
to an equilibrium current, known as persistent currenti 

Observations of persistent currents in metallic and 
semiconducting rings^ stimulated many theoretical pa- 
pers focused on the electron-electron (e-e) interaction and 
disorder in such systems^. A quantitative understanding 
of the observed amplitude of the currents is so far not sat- 
isfactory and the role of the e-e interaction and disorder 
in mesoscopic rings need further researchi^. In this work 
we focus on the persistent current of interacting spinless 
electrons in a ID ring with a single scatterer. Let us 
review recent results for this simplified problem. 

If the e-e interaction is ignored, the leading behavior 
of the persistent current (J) as a function of the mag- 
netic flux (0) and ring circumference (L) can be derived 
analytically for an arbitrary scatterer—. It is a periodic 
function of (j) with period 4>o — h/e, it decays as 1/L, 
and its shape is determined by the absolute value of 
the transmission amplitude tk F of the scatterer at the 
Fermi wave vector hp. In the limit of small tk F it reads 
I cx L~ 1 \tk F | sin(2w(f>/(f>o). These results were derived 
for a continuous ring£ and hold also for the ring-shaped 
tight-binding lattice at half filling 6 . 

For a repulsive e-e interaction in the continuous ring, 
the spinless persistent current can still be derived ana- 
lytically in the Luttinger liquid modem For L — > oo and 
for an arbitrary strength of the scatterer, the current is 
/ cx i - " -1 sin(27r</>/0o), where the power a > is uni- 
versal, depending only on the e-e interaction. 

The authors of Ref. @ obtained the spinless persistent 
current microscopically by solving the Hubbard chain 



with nearest-neighbor hopping and interaction. Using 
the renormalization group, they indeed found for a ID 
ring with a single scatterer the current which decays 
faster than 1 /L and is a sine- like function of 0. The shape 
/ cx sin(27r</>/(/>o) and the universal scaling / cx L~ Q_1 
were confirmed for long chains and/or strong scatterers. 

They also concluded that the sine-like shape and the 
scaling can only be obtained in the correlated 
many-body model, not in the Hartree-Fock approxima- 
tion ignoring the correlation effects. They used a non- 
self-consistent Hartree-Fock approximation and, indeed, 
found the qualitatively wrong results. The self-consistent 
Hartree-Fock approximation, according to these authors, 
drives the model with a single scatterer into the charge- 
density-wave groundstate with a finite single-particle 
gap, which is wrong because a single scatterer cannot 
change the bulk properties of the system. Hence, no 
attempt has been done to test the dependence / cx 
£-a-i sm (27r0/0o) in the self-consistent Hartree-Fock 
approximation. However, there is a motivation to do so. 

Indeed, recently the self-consistent Hartree-Fock ap- 
proximation has been used to study tunneling of the 
weakly-interacting electron gas through a single scat- 
terer in a continuous ID wire with leads. It has been 
found that there is no charge-density-wave groundstate 
in such model, moreover, the tunneling probability has 
been found to decay with the wire length as L~ 2a , 
in accord with correlated models 8,9 . All this suggests 
that an attempt to obtain the persistent current I cx 
L _a_1 sin(27r0/0o) m the self-consistent Hartree-Fock 
approximation might be meaningful. 

We present such attempt in our present work. Using 
the self-consistent Hartree-Fock approximation at zero 
temperature, we study the persistent current of the in- 
teracting spinless electron gas in a continuous ID ring 
with a single S barrier. We compare our results with the 
results of correlated models, with the Luttinger liquid 
models and the microscopic Hubbard model 6 -. A good 
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agreement is found. First, the persistent current decays 
faster than L~ x and eventually like L~ Q_1 , where a > 
is universal. Second, the persistent current is a sine-like 
function of magnetic flux. 

The paper is organized as follows. In Sect. II we dis- 
cuss the Hartree-Fock model of the ID ring. In Sect. Ill 
we present the results for the Hartree-Fock potentials. In 
Sect. IV the Hartree-Fock results for persistent currents 
are presented and compared with the results of correlated 
models. A summary is given in Sect. V. 



where fcj is the quantum number of the i-th electron. The 
wave functions tpk {%) obey the Hartree-Fock equation 



2m 



d 2tt 



dx L cj)Q 
+ U H {x) + U F (k,x) 
with the boundary condition 



7<5(x) 



ipk(x) = e k ip k {x) (6) 



II. HARTREE-FOCK MODEL OF THE ID RING 

We consider TV interacting ID electrons with free mo- 
tion along a circular ring threaded by magnetic flux 
4> = BS = AL, where S is the area of the ring, B is 
the magnetic field (constant and perpendicular to the 
ring area), and A is the magnitude of the resulting vec- 
tor potential (circulating along the ring circumference). 
Such electron system is described by Hamiltonian 



H = 



N 

E 



h 2 ( . d 2tt r . 



2m 



3 L 

N 



dxj 



(1) 



where m is the electron effective mass, xj is the coordi- 
nate of the j-th electron along the ring, V(xj — xt) is the 
e-e interaction, and "f5(x) is the potential barrier due to 
the scatterer. The eigenfunction ^{x\,X2, ■ ■ ■ ,xn) and 
eigenenergy E obey the Schodinger equation 



tp k (x + L) = ip k {x) 



(7) 



where 



U H (x) =J2j dx'V{x-x')\^{x')\ 2 (8) 



is the Hartree potential and 



U F (k, x) 



f dx'V{x-x')^ k {x')rw{x')^{x) (9) 



is the Fock nonlocal exchange term (expressed as an effec- 
tive potencial for further convenience). In the equations 
(JSJl and 10 we sum over all occupied states k'. 
We introduce the wave functions ip k (x) by substitution 



* 7,(.r) = y /,(•'■) exp ( -i^J-X 



(10) 



Equations JBJ-J^J then give the Hartree-Fock equation 



= Et> 

with the cyclic boundary condition 



(2) 



^(xi, . . . ,X{ + L, . . .,x N ) = *(^!, . ..,Xi, .. -,x N ) (3) 

for i = 1, 2, . . . , N. In the groundstate with the eigenen- 
ergy Eq the persistent current reads 



(4) 



In the Hartree-Fock approximation, the wave function 
^ is approximated by the Slater determinant 



. . -,x N ) 



^ki{xx) 



Tpk x {xN) 



ip kN (xi) ■■■ il> kN (x N ) 



, (5) 



1LJ*L 

2m dx 2 



+ j5(x) + U H {x) + U F {k, x) 



with the boundary condition 



Vk{x) 
£Wk(x) (11) 

x) , (12) 



where the potentials Uh(x) and U F (k,x) are still given 
by equations JHJ and ©, but with -0 replaced by ip. 

From equations ifTjl- lfP^I) the groundstate energy E — 
(^\H\^>) can be expressed after simple algebra as 



(fk(x + L) = exp ( «27r— I tp k ( 



£fe - ^ (fk \Uh{x) + U F (k,x)\ ipk) 



(13) 



Setting this expression into Q we obtain the persistent 
current in the Hartree-Fock approximation. Once E^{<j)) 
is known, dE$(<j)) / d<f> can be calculated numerically. 
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As in Refs. ^] and fill we simplify the equation JSJ) as 
U F (x) ~ - [ dx'V(x - x')Re (z')Vv (a;)] (14) 



by noticing that J2k' ^k'( x ')^k'(x) — 5(x — x'). This 
'almost closure relation' was tested in Fig. 2 of Ref. ITU 
so we do not repeat the same test here. Our comments on 
the reliability of approximation 114|l are given at the end 
of Sect. IV. Unlike the exact form ©, the Fock potential 
(|14|l is local and independent on k. This saves computer 
time and memory and allows to study long rings. 

In absence of the S barrier, the solution of equation 
© is i^k( x ) = e%kx an d the potentials Uh{x) and Up(x) 
become ^-independent constants 7 -. These constants, Uh 
and Uf, only introduce a rigid shift of the energy scale 
which has no physical effect. Therefore, in the following 
text we simply consider the potentials Uh{x) and Uf(x) 
reduced by constants Uh and Uf, respectively. 

We now discuss the solution of equations and 
(|12(l assuming that the potentials Uh(x) and Uf{x) are 
known. Consider first the ring region x E (— L/2, L/2) as 
a straight-line segment of an infinite ID wire. Inside the 
segment the potential is r y5(x) + Uh{x) + Up(x), outside 
we keep it zero. Therefore, the wave function outside is 



be- lkx , x<-L/2, 



Vk (x) = ce tkx + de- lkx , x > L/2 . 
The amplitudes a and b are related to c and d by 



To 



(15) 



(16) 



(17) 



where Tn is the transfer matrix 



T 



j_ 

'I 

-Lk 



(18) 



with tk and r k being the transmission and reflection am- 
plitudes of the electron impinging the region (—L/2, L/2) 
from the left. To come back to the ring threaded by mag- 
netic flux, we relate (fk(—L/2) and ip k (L/2) (eqs. IT5l[T^|) 
through the boundary condition 1|12|) . Combining this 
relation with eqs. I|17|l and (|18f) we obtain the equation 



where 



T 



T 



a \ , q> , I a 

exp(i27r — ) 

4>o 



J_„ikL r k r ikL 



-ikL 



(19) 



(20) 



Thus exp(i2TT<p/<po) is the eigenvalue of matrix T. The 
product of the eigenvalues of a matrix is given by its de- 
terminant, which is unity in this case. The second eigen- 
value is thus exp(— «27T</)/(/)o). Their sum is equal to the 
matrix trace, which gives the equation for the spectrum, 5 



cos 



H) 


= Re 


exp(— ikL) 






t k 



(21) 



Numerical solution of equation (|21|l has to be combined 
with computation of the transmission amplitude tk ■ We 
compute tk and rk in the same way as described for the 
open wire 7 . The solution of equation (|21|l gives us the 
dependence k(<f>) and finally Sk(4>) — ft 2 k 2 ((ft) /2m. 

For each k(cj>) we calculate the wave function ifik(x) as 
follows. From eqs. I|19|) an< 4 l|2U|) we obtain the amplitude 



— - — e 

rk r k 



i{2ir<j>/4, +kL) 



(22) 



and the amplitude b by normalizing the wave function. 
Then we express from equation (|15f) the boundary condi- 
tions ifk(—L/2) and dfk(~L/2)/dx. We solve the equa- 
tion (flip with these boundary conditions by using the 
same finite-difference method as described for open wire 7 . 

Once we obtain (fik(x), we can also perform the self- 
consistent Hartree-Fock calculation. In the first iteration 
step we solve equation (|ll|l for the non- interacting gas. 
The resulting (fik(x) is used to evaluate the Hartree and 
Fock potentials. These potentials are used in the second 
iteration step to obtain new <fk(x) and new potentials 12 , 
etc., until the energies £k do not change any moreii. 

In this work we apply the above Hartree-Fock scheme 
to the GaAs ring with electron density n = 5 x 10 7 m — 1 , 
effective mass m = 0.067 mo, and e-e interaction 



V(x-x') = V e-\ x - x '\ /d , 



(23) 



where Vq = 34 meV and d = 3 nm. We adopt the finite- 
ranged interaction l|23|l due to the following reasons. 
The e-e interaction in electron gas is always screened, 
so a reasonably chosen screened e-e interaction improves 
the Hartree-Fock approximationii The e-e interaction 
(12311 has already been used to study the many-body ID 
systems 14 , it has been showni to mimic the screened e-e 
interaction in a realistic ID system quite well. Finally, 
we want to compare our study with the correlated models 
which also rely on the e-e interaction of finite range. 



III. HARTREE-FOCK POTENTIAL 

In the following we label the transmission and reflec- 
tion amplitudes of the bare 5 barrier as tk and fk- It 
holds thai*! 5 . t k = k/(k + i() and r k = -i(/(k + <), 
where C, = jm/h 2 . Since kp and m are fixed, instead of 
using 7 we parametrize the 5 barrier by its transmission 

I ~ 1 2 

coefficient \t kF ■ 
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FIG. 1: Self-consistent Hartree-Fock potential Uh{x) + Uf{x) 
along a ID ring, induced by the (not shown) barrier 'yS(x) due 
to the scatterer at x — 0. The 8 barrier is adjusted to have the 
bare transmission \tk F \ = 0.851 at the Fermi level (14 meV). 
The length of the ring, L, is 20 /im (in the figure only a small 
region around the scatterer is shown). The solid curve in fact 
involves two overlapping curves for the rings with magnetic 
flux 4> — * + and <j> = O.250o- Notice that we find the same 
curve for any cj> as long as we are not far from the 8 barrier (see 
the next figure). For comparison, the dashed curve shows the 
result of Ref. the self-consistent Hartree-Fock potential 
in a GaAs wire connected to leads, with the 8 barrier and 
parameters as in our ring. 



The solid curve in figure ^ shows the self-consistent 
Hartree-Fock potential in the ring, induced by the scat- 
terer at x — 0. The ring is 20 long, we show only a 
small region around the scatterer. The solid curve in fact 
involves two curves overlapping each other, one for the 
ring with very small flux <fi — > + and one for the ring with 
4> = O.250O- The dashed curve is the result taken from 
figure 2 of Ref. 0. It is the self-consistent Hartree-Fock 
potential in the ID GaAs wire with the same scatterer 
and parameters as in our ring. All curves exhibit the 
Friedel oscillations with period Af/2. Since we compare 
the long ring and long wire in a close vicinity of a; = 0, 
the agreement of the solid and dashed curves is not sur- 
prising. As shown in the next figure, all three curves 
differ remarkably at large distances from the scatterer. 

Figure [3 shows all three curves along the entire half 
length. Here the Friedel oscillations are too dense to 
be distinguishable, but we can observe their asymptotic 
decay. We see in the top and middle panel, that the 
amplitude of the Friedel oscillations in the ring saturates 
at large distances from the scatterer. This saturation is 
due to the circular shape and was already observed in 
the pioneer Hartree-Fock study of the continuous ring 10 . 

However, the authors of Ref. discussed also the wire 
with leads, but without performing any self-consistent 
Hartree-Fock study. They estimated analytically for ex- 
tremely strong e-e interaction, that in the infinite wire 
the amplitude of the Friedel oscillations saturates with 
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FIG. 2: The same Hartree-Fock potentials as in the preceding 
figure, but along the entire half length. The top panel and 
middle panel show the results for the ring with <f> — > + and 
4> = O.250o, respectively. For comparison, the bottom panel 
shows the result for the GaAs wire with leads, taken from 
Ref. Inset to each panel shows the same data in log scale, 
separately for the positive (top inset) and negative (bottom 
inset) data. Inset to the bottom panel shows, that the decay of 
Friedel oscillations in the wire is linear in log scale. It follows 
the straight (dashed) line with slope /3 = 0.9. A detailed 
analysis 7 shows that j3 —* 1 at extremely large distances. 



distance similarly as in the ring. In other words, at 
large distances from the scatterer Uh(x) + Uf(x) cx 
x~P cos(2kpx + const), but with (3 — 0. This is the 
charge-density-wave groundstate. In Ref. this ground- 
state was concluded to hold for any e-e interaction 
strength, but without providing explicit self-consistent 
Hartree-Fock results. Nevertheless, due to this charge- 
density-wave groundstate, the self-consistent Hartree- 
Fock method was considered to be wrong 6 . 
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Recently 7 ., the self-consistent Hartree-Fock study has 
been applied to the continuous ID wire with a single scat- 
terer, connected to leads. It has been shown, that the 
charge-density-wave groundstate does not exist in that 
model at least for weak e-e interaction. In the bottom 
panel of figure we show the result of Ref. in order 
to compare with the Hartree-Fock potential in the ring. 
Inset to the bottom panel shows, that the asymptotic de- 
cay of Friedel oscillations in the wire is linear in log scale. 
This implies the dependence oc with positive j3 > 0. 
A detailed analysis 7 shows that (3 — > 1 in the limit of very 
long wires, i.e., there is no tendency to a non-decaying 
((3 = 0) periodic wave. Saturation of the oscillation am- 
plitude in the ring is due to the ring geometry and does 
not mean any charge-density-wave groundstate. 



IV. RESULTS FOR PERSISTENT CURRENT 

For noninteracting spinless electrons, the persistent 
current in the continuous ID ring with a single scatterer 
can be derived analytically 5 . For even N (considered 
throughout this paper) the result reads 

^e^arc cosd^lcos^-Tr ]) , n 

where 4>' = 2n<t>/(j)Q. Equation l|24|l holds also for the 
tight-binding model at half-filling^. For small tk F 



I=^\t kF | sinful 



2L 



(25) 



For the repulsive e-e interaction, as already mentioned, 

J oc L'"- 1 sin(0') (26) 

within the Luttinger-liquid model—. Due to the e-e inter- 
action the current vanishes faster than l/L. 

These analytical predictions were tested in a numer- 
ical "experiment" of Ref. ;6. The authors of that 
work studied a ring-shaped Hubbard lattice with nearest- 
neighbor hopping and nearest-neighbor e-e interaction. 
They found for the ID ring with a single scatterer, that 
the current indeed vanishes faster than 1 / L and is a sine- 
like function of <fr' . Their calculation was based on the 
renormalization-group (RG) technique allowing a highly 
reliable microscopic treatment of correlations. 

Let us compare the RG calculation of Ref. |6| with our 
Hartree-Fock model. In panel a of figure [3| we show the 
persistent current LI((f>') calculated by our self-consistent 
Hartree-Fock approach, in panel b we show the RG re- 
sult of Ref. 6 in terms of NI(<j)'). For the purpose of 
comparison we follow Ref. |6| and normalize our LI((f>') 
data by cvf/2. Our ring is continuous and the ring of 
Ref. la is the Hubbard lattice at half-filling. Therefore, 
we cannot expect quantitative agreement. However, we 
see excellent qualitative agreement. In both cases the e-e 
interaction preserves the sine-like dependence on cj)' and 
in both cases the current decays with L faster than l/L. 
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FIG. 3: Panel a shows the persistent current LI(4>), calculated 
by our self-consistent Hartree-Fock approach for \tk F = 0.64 
and various L. For qualitative comparison, panels b and c 
show the data taken from Fig. 3 of Ref. [(J These data were 
obtained— for a ID ring-shaped Hubbard lattice with the same 

I ~ 1 2 

\tk F and same N as in panel a, where N = 32, 64, . . . for 
L — 0.64/im, 1.28/im, . . . , respectively. In particular, panel 
b shows the results of the RG solution while panel c shows 
the results of the (non-self-consistent) Hartree-Fock solution. 
The RG solution accounts for the effect of correlations. 



The authors of Ref. G also performed their own (non- 
self-consistent) Hartree-Fock calculation. For the pur- 
pose of comparison, the results of that calculation are 
shown in panel c of Fig. |3J One can see that these re- 
sults exhibit non-physical deviations from the sine-like 
shape. This non-physical behavior was ascribed- to the 
absence of correlations in the Hartree-Fock model, but 
we can now ascribe it to the absence of self-consistency 
(our self-consistent Hartree-Fock results in panel a clearly 
exhibit the sine- like shape). 
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FIG. 4: Panel a shows the persistent current LI (ft), calcu- 
lated by our self-consistent Hartree-Fock approach for vari- 

I ~ i 2 

ous L and \tk F \ = 0.22. For qualitative comparison, panel b 
shows the RG data taken from Fig. 4 of Ref. [f| These data 
were obtained— for a ID ring-shaped Hubbard lattice with the 

I ~ 1 2 

same \tt F \ and same N as in panel a, where N = 32, 64, . . . 
for L = 0.64pm, 1.28pim, . . . , respectively. 
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In figure 01 we again compare our Hartree-Fock results 
with the RG results of Ref. 0, but for a stronger S barrier. 
We see a similar qualitative agreement as in Fig. |3 the 
peak of the current is now shifted closer to 4>' = n/2. 

In figure [5] we present the persistent current LI(<f>') 

I ~ i 2 

for the 6 barrier as strong as \tk F \ = 0.03. Our self- 
consistent Hartree-Fock results are shown in panel a. We 
want to compare these results with the scaling law due 
to the Luttinger liquid model (eq. (J2SJ)). To evaluate 
the scaling law Q26f) . we reformulate it as follows^. The 
idea is to replace the bare transmission amplitude t^ F in 
the non-interacting scaling law l|25() by the transmission 
amplitude of the interacting electron gas 9 , 



FIG. 5: Persistent current LI(4>') for tfep = 0.03 and var- 
ious L. Panel a shows our self-consistent Hartree-Fock data, 
panel b the data typical of the Luttinger liquid (see the text). 



the right hand side of (|27|l holds for small tk F and/or 
large L. For small tu F one indeed obtains the scaling law 



/ = luL- 



sm 



(28) 



where w = evF \tk F \ d a /2 \fk F \- The formula l|27|l was 
derived^ by means of the RG method in the limit of weak 
e-e interaction (a « 1). In that limit a is given by^ 



t kF (d/LY 



r kF \ 2 + \t kF \ 2 (d/L)2« 



tk F 
\fk F 



(d/Lf, (27) 



where d is the range of the e-e interaction V(x — x') and 



V(0) - V(2k F ) 



(29) 



where V(q) is the Fourier transform of the e-e interac- 
tion V(x — x'). The Fourier transform of our interaction 




L [urn] 

FIG. 6: Persistent current LI(cj>' = vr/2) as a function of L 
for various S barriers. Our Hartree-Fock results are shown by 
open symbols connected by dashed lines. The full lines in the 
left panel show the scaling law 12811 with power a — 0.0855. 
In this calculation the prefactor ui in eq. 1281 was adjusted to 
fit the Hartree-Fock data quantitatively. The full lines in the 
right panel show the exact dependence 1241 with tk F replaced 
by the exact amplitude (1271 . 

(|23J) reads V(q) = 2V d/(l + q 2 d 2 ). Wc set this for- 
mula and our parameters into equation l|29[l . We obtain 
a = 0.0855. We evaluate the formula (|28[l and we show 
this result in panel b of figure [5J Clearly, this result and 
the Hartree-Fock result of panel a are in excellent quali- 
tative and good quantitative agreement. 

In figureEJwe show the persistent current LI(cf)' = it/2) 
as a function of L for various S barriers. Our Hartree- 
Fock results are shown by open symbols connected by 
dashed lines. The full lines in the left panel show the 
scaling law (|2*S1) with the power a — 0.0855 obtained 
from the universal formula l|29[) and with the factor lu 
adjusted to fit the open symbols quantitatively. We see 
that the Hartree-Fock data reproduce the universal law 
/ oc L~ a ~ l in the limit of small tk F - Of course, the power 
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FIG. 7: Persistent current LI((f>' = tt/2) versus L for various S 
barriers and various e-e interaction strength (various a) . Our 
Hartree-Fock results are shown by open symbols connected by 
dashed lines. The full lines show the scaling law 1281 with the 
factor ui adjusted to fit the Hartree-Fock data quantitatively. 
In this calculation the powers a = 0.0277 and a — 0.0561 
correspond to the the e-e interaction amplitudes Vo = 11 meV 
and Vo = 22.3 meV, respectively. Note that we can keep the 
same value of a by using the e-e interaction 123H with many 
different combinations of Vo and d. However, as long as a 
remains the same, the choice of Vo and d has no effect on our 
results (Hartree-Fock potentials, eigenenergies, currents). 



law l|28|) is asymptotic and cannot agree with the Hartree- 
Fock data for any tk F ■ A good agreement is obtained if 
we plot directly the exact formula l|24|) with tk F replaced 
by the exact amplitude (|27|l . This exact dependence is 
shown in a full line in the right panel of figure |SJ The full 
lines and Hartree-Fock data now agree quite well for any 
tk F and scale as / oc i""^ 1 for small tk F . Why a small 
quantitative difference appears is explained below. 

Eventually, for small enough tk F our Hartree-Fock data 
always follow the scaling law / oc L~ a ~ x . This is demon- 
strated in Fig[7|for several tk F and a. 

Now we add a few arguments on the reliability of ap- 
proximation (|14f> . Our Hartree-Fock results systemat- 
ically agree with the results of correlated models. This 
would hardly be possible in a wrong Hartree-Fock model. 
In particular, we have recovered the power law / oc 
L - "^ 1 , where a is given by universal expression 12911 . 
It is highly improbable that such nontrivial physical law 
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could be a purely accidental consequence of wrong ap- 
proximation in the Hartree-Fock model. 

Since it relies^ on the 'almost closure relation' 
J2k< Tpk' ( x ')V'fc' ( x ) — 8{x—x'), approximation (|14|l in fact 
underestimates the exact Fock potential 10 and overes- 
timates the electron transmission and persistent current. 
Using of the exact Fock potential J^J would just remove 
a small quantitative difference between the Hartree-Fock 
and correlated results in the right panel of Fig. HJ 

Why the approximation (|14|l works so well can be qual- 
itatively understood as follows. It can be seen to be a 
perfect approximation of equation @ in the extreme case 
V(x — x') oc 5(x — x'). Thus, the shorter the range of the 
interaction V(x — x') the better the approximation (|14|l . 
Our interaction l|23|) is short-ranged and in this respect 
more appropriate than the bare Coulomb interaction for 
which the approximation (|14|l seems to work wellfi^iii 
For V(x— x') oc 5(x — x') the Hartree and Fock potentials 
cancel to zero, the exponentially decaying e-e interaction 
(|23[l still produces a weak Hartree-Fock potential. 

V. CONCLUDING REMARKS 

In conclusion, using the self-consistent Hartree-Fock 
approximation at zero temperature, we have studied the 
persistent current of the weakly-interacting spinlcss elec- 
tron gas in a continuous ID ring with a single 8 barrier. 
We have compared our results with the results of corre- 



lated models, with the Luttinger liquid models and the 
microscopic Hubbard model&. We have found a good 
agreement. In particular, in the limit of strong 8 bar- 
riers our self-consistent Hartree-Fock approximation re- 
produces the scaling of the persistent current in the form 
/ oc L _Q_1 sin(27r0/0o), where a is a universal power, 
depending only on the e-e interaction. We can therefore 
conclude, that the universal power law is not exclusively 
the consequence of correlations; it can also be obtained 
in the self-consistent Hartree-Fock approximation. 

We believe that the self-consistent Hartree-Fock ap- 
proximation would confirm the L~ a ~ l scaling also in the 
limit of weak 8 barriers, where the L _Q_1 scaling arises 
only at extremely large L. Unfortunately, our present 
computer equipment is too slow for this task. 

The self-consistent Hartree-Fock approximation has 
recently! been used to study tunneling of the weakly- 
interacting electron gas through a single scatterer in a 
continuous ID wire with leads. It has been foundi in the 
limit of strong 8 barriers, that the tunneling probabil- 
ity decays with the wire length as L~ 2a , in accord with 
correlated models^. A quite similar agreement with cor- 
related models shows our present Hartree-Fock study by 
reproducing the L" Q_1 scaling of persistent currents^. 
Acknowledgments 

M.M. was supported by the APVT grant APVT-51- 
021602, R.N. by the VEGA grant 2/3118/23. 



* Electronic address: martin.mosko@savba.sk 

1 Y. Imry, Introduction to Mesoscopic Physics (Oxford Uni- 
versity Press, Oxford, UK, 2002). 

2 L. P. Levy, G. Dolan, J. Dunsmuir, and H. Bouchiat, Phys. 
Rev. Lett. 64, 2074 (1990); V. Chandrasekar, R. A. Webb, 
M. J. Brady, M. B. Ketchen, W. J. Gallagher, and A. 
Kleinsasser, Phys. Rev. Lett. 67, 3578 (1991); D. Mailly. 
C. Chapelier, and A. Benoit, Phys. Rev. Lett. 70, 2020 

(1993) ; E. M. Q. Jariwala, P. Mohanty, M. B. Ketchen, 
and R. A. Webb, Phys. Rev. Lett. 86, 1594 (2001); W. 
Rabaud, L. Saminadayar, D. Mailly, K. Hasselbach, A. 
Benoit, and B. Etienne, Phys. Rev. Lett. 86, 3124 (2001). 

3 Reviews of this issue can be found in S. Viefers, P. Koski- 
nen, P. S. Deo, and M. Manninen, Physica E 21, 1 (2004), 
in U. Eckern, and P. Schwab, J. Low Temp. Phys. 126, 
1291 (2002), and in Ref.fl 

4 Y. Imry and R. Landauer, Rev. Mod. Phys 71, S306 
(1999). 

5 A. O. Gogolin and N. V. Prokof'ev, Phys. Rev. B 50, 4921 

(1994) . 

6 V. Meden and U. Schollwock, Phys. Rev. B 67, 035106 
(2003). 

7 M. Mosko, A. Gendiar, P. Vagner, and Th. Schappers, sub- 



mitted to Phys. Rev. B. 

8 C. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 68, 1220 
(1992). 

9 K. A. Matveev, D. Yue, and L. I. Glazman, Phys. Rev. 
Lett. 71, 3351 (1993); D. Yue, L. I. Glazman, and K. A. 
Matveev, Phys. Rev. B 49, 1966 (1994). 

A. Cohen, K. Richter, and R. Berkovits, Phys. Rev. B 57, 
6223 (1998). 

1 A. Cohen, R. Berkovits, and A. Heinrich, Int. J. Mod. 
Phys. B 11, 1845 (1997). 

2 Inspired by the Schrodinger/Poisson solver for inversion Si 
layers [T. Ando, A. B. Fowler, and F. Stern, Rev. Mod. 
Phys. 54, 437 (1982)], we use a similar iterative trick. 

3 N. W. Ashcroft, and N. D. Mermin,Sohd state physics, 
Sounders College Publishing, Ed. D. G. Crane, Cornell 
University, USA 1976. 

4 V. Meden, W. Metzner, U. Schollwock, O. Schneider, T. 
Stauber, and K. Schonhammer, Europhys. J. B 16, 631 
(2000). 

J. H. Davies, The Physics of Low-Dimensional Semicon- 
ductors: An introduction (Cambridge University Press, 
Cambridge, UK, 1998). 



